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Abstract. It has recently been shown that f(R) theories formulated in the Palatini variational 
formalism are able to avoid the big bang singularity yielding instead a bouncing solution. The 
mechanism responsible for this behavior is similar to that observed in the effective dynamics 
of loop quantum cosmology and an f(R) theory exactly reproducing that dynamics has been 
found. I will show here that considering more general actions, with quadratic contributions of 
the Ricci tensor, results in a much richer phenomenology that yields bouncing solutions even in 
anisotropic (Bianchi I) scenarios. Some implications of these results are discussed. 



1. Introduction. 

In the last years it has been shown that though many f{R) models have the ability to produce 
late-time cosmic acceleration and fit well the background expansion history, they are not in 
quantitative agreement with the structure and evolution of cosmic inhomogeneities (see and 
references therein). Additionally, the fact that matter is concentrated in discrete structures like 
atoms causes the modified dynamics to manifest also in laboratory experiments [21 [3l which 
confirms earlier suspicions on the viability of such models according to their corresponding 
Newtonian and post-Newtonian properties [5]. This is a very disturbing aspect of the models 
with infrared corrections, which demands the consideration of a microscopic description of the 
sources and prevents the use of macroscopic, averaged representations of the matter. A careful 
analysis of this point put forward the existence of non-perturbative effects induced by the Palatini 
dynamics in a number of contexts [Q [21 [71 13| . Though this certainly is an undesired property 
of infrared-corrected models, it could become a very useful tool for models with corrections 
at high curvatures. Can we construct singularity- free cosmological models that recover GR at 
low curvatures using the non-perturbative properties of Palatini theories? As we report here, 
ultraviolet-corrected Palatini models turn out to be very efficient at removing the big bang 
cosmic singularity in various situations of interest. 



2. Non-singular f{R) cosmologies. 

The recent interest in the dynamics of the early-universe in Palatini theories arose, in part, 
from the observation that the effective equations of loop quantum cosmology [8] (LQC) could 
be exactly reproduced by a Palatini f{R) Lagrangian [9]. In LQC, non-perturbative quantum 
gravity effects lead to the resolution of the big bang singularity by a quantum bounce without 
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introducing any new degrees of freedom. Though fundamentahy discrete, the theory admits a 
continuum description in terms of an effective Hamiltonian that in the case of a homogeneous 
and isotropic universe filled with a massless scalar field leads to the following Friedmann equation 
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where pcnt ~ ^-^^PPianck- At low densities, p/pcrit ^ the background dynamics is the same 
as in GR, whereas at densities of order pcrit the non-linear new matter contribution forces the 
vanishing of and hence a cosmic bounce. This singularity avoidance seems to be a generic 
feature of loop-quantized universes [10] . 

Palatini f{R) theories share with LQC the fact that the modified dynamics that they produce 
is not due to the existence of new dynamical degrees of freedom but rather to non-linear effects 
produced by the matter sources, which contrasts with other approaches to quantum gravity and 
to modified gravity. This allows to relate Eq.([TJ with the corresponding f{R) equation 
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Taking into account the trace of the field equations (see p] for details), which for a massless 
scalar becomes RJb. — 2/ = 2k^p and implies that p = p{R), one finds that a Palatini f{R) 
theory able to reproduce the LQC dynamics ([T|) must satisfy the differential equation 
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where A = ^2{RfR - 2f){2R, - [R/r - 2/]), B = 2^R,fR{2RfR - 3/), and R, = k^c- 
Imposing suitable boundary conditions, one can show [£] that the function 
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faithfully captures the LQC dynamics from low curvatures up to the non-perturbative bounc- 
ing region. Though the function Q implies that the LQC Lagrangian is an infinite series, 
which is a manifestation of the non-local properties of the quantum geometry, one can find 
non-singular cosmologies with a finite number of terms. In fact, a simple quadratic Lagrangian 
of the form f{R) = R + R? / Rp does exhibit non-singular solutions for certain equations of state 
[IT] depending on the sign of Rp. To be precise, if Rp > the bounce occurs for sources with 
w = P/ p > 1/3. If Rp < 0, then the bouncing condition is satisfied by w < 1/3. 



Besides avoiding the development of curvature singularities, bouncing cosmologies can solve 
the horizon problem |I2], which makes them interesting as a substitute for infiation. However, 
it has been found |13j that f{R) models that develop a bounce when the condition //j = is 
met turn out to exhibit singular behavior of inhomogeneous perturbations in a flat, dust-filled 
universe. Further insight on the robustness of the bounce under perturbations was obtained 
in [13] studying the properties of f{R) theories in anisotropic spacetimes of Bianchi-I type, for 
which ds^ = -dt^ + E?=i aj{t){dx' f. In space-times of this type with a perfect fluid, one can 
derive a number of useful analytical expressions. In particular, one finds that the expansion 

6 = Hi and the shear = Yl,i (^Hi — |^ are given by 



where the constants Cij = —Cji set the amount and distribution of anisotropy and satisfy the 
constraint C12 + C23 + C31 = 0. In the isotropic case, Cij = 0, one has cr^ = and 9^ = 9H^, 
with given by Eq.([2]). Now, since homogeneous and isotropic bouncing universes always 
require the condition //^ = at the bounce |l4j, a glance at ([5]) indicates that the shear diverges 
as ~ This shows that any isotropic f{R) bouncing model will develop divergences when 

anisotropies are present. Moreover, one can check by direct calculation that the Kretschman 
scalar Rf.uapR'""'''' = KEii^i + Hff + EfH^ + HfHl + H^Hl) diverges at least as ~ 
which is a clear geometrical pathology and signals the presence of a physical singularity. 



3. Nonsingular cosmologies beyond f{R). 

The previous section provides reasons to believe that Palatini f{R) models are not able to 
produce a fully satisfactory and singularity-free alternative to GR. Though the homogeneous 
and isotropic case greatly improves the situation with respect to GR, the existence of divergences 
when anisotropies and inhomogeneities are present spoils the hopes deposited on this kind of 
Lagrangians. To the light of these results, we explored new Palatini theories [II] to determine 
if the introduction of new geometrical invariants in the gravitational action could avoid the 
problems that appear in the f{R) models. We considered Lagrangians of the form f{R, Q), with 
Q = Rfj^ijR^'^ . Using the particular Lagrangian 
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it was found that completely regular bouncing solutions exist for both isotropic and anisotropic 
homogeneous cosmologies filled with a perfect fluid. In particular, we found that for a < the 
interval 0<t(;<l/3is always included in the family of bouncing solutions, which contains the 
dust and radiation cases. For a > 0, the fluids yielding a non-singular evolution are restricted 
which implies that the radiation case w = 1/3 is always nonsingular. For a 
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detailed discussion and classification of the non-singular solutions depending on the value of the 
parameter a and the equation of state w see [11] . 

The field equations that follow from the Lagrangian ([6]) when R^^j is assumed symmetric were 
derived in 1151 and take the form 
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where /r = duf and fq = dqf . The connection equation in ([7]) can be solved introducing an 
auxiliary metric haji such that it takes the form \/—hh'^'^ = 0, which implies that F^^^ can 
be written as the Levi-Civita connection of hi^^. When the matter sources are represented by a 
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one can show that h^y and h^^^ are given by |15j 
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where \ = ^ ^^P + I + ^, A, = ^2^A + ^, A2 = [x ± ^X'-k\p + P) 

= [Ai(Ai — A2)]^''^. The metric field equation in ([7]) takes the following form 
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In this expression, the functions /, Ai, and A2 are functions of the density p and pressure 
P. In particular, in a universe filled with radiation, for which R = 0, the function Q in our 
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quadratic model boils down to Q = 



This expression recovers 
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the GR value at low curvatures, Q ~ 4:{k^p)'^/3 + 32(«;^p)^/9i?p + . . . but reaches a maximum 
Qmax = 3i?p/16 at K^Pmax = 3i?p/16. At /9max the shear also takes its maximum allowed value, 

namely, cr'^ax — \/3/16i2p ^(Cx2 + C23 + C^i), which is always finite, and the expansion vanishes 
producing a cosmic bounce regardless of the amount of anisotropy. 

It should be noted that the choice of a symmetric Ricci tensor in the analysis of f{R, Q) bounc- 
ing cosmologies presented above is not arbitrary. As pointed out in [16j, the antisymmetric 
part of the Ricci tensor introduces new dynamical degrees of freedom in the form of a massive 
vector field (see also [T7j for a related result). If one looks for a framework suitable for the 
description of the effective dynamics of LQC (including anisotropics) and, more generally, of 
other theories of quantum geometry not involving new degrees of freedom, it seems natural to 
impose constraints on the spectrum of possible Lagrangians to avoid new propagating fields. In 
this sense, it should be noted that the f{R, Q) theories discussed here are also able to reproduce 
[18j other aspects of the expected phenomenology of quantum gravity at the Planck scale. In 
particular, without imposing any a priori phenomenological structure, the quadratic Palatini 
model ([6]) predicts an energy-density dependence of the metric components that closely matches 
the structure conjectured in models of Doubly (or Deformed) Special Relativity and Rainbow 
Gravity |19] . This confirms that Palatini theories represent a new and powerful framework to 
address different aspects of quantum gravity phenomenology. 
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